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Abstract 

An alternative model of the fragmentation of the Lund string, investi- 
gating the hypothesis of helix-like ordered gluon field, is compared with 
hadronic Z° data. A significant improvement in the description of various 
measured quantities is achieved. In particular, the existence of correlations 
between longitudinal and transverse components of hadron momenta (as sug- 
gested by the helix string model) seems to be supported by the data. 

1 Introduction 

The possibility that a helix-like ordered gluon field emerges at the end of a parton 
cascade was first discussed in JT), on the basis of a study of optimal packing of 
soft gluons in the phase space. In their pioneering paper, the authors suggested to 
study a helix parametrization related to the rapidity difference along string. The 
search of such a helix string signature was performed but no significant signal was 
found 0. 

In this paper, we will study an alternative parametrization of the helix string, 
which may be at the origin of some persistant discrepancies between data and mod- 
elling. 

The paper is organized as follows : in section 2, a very brief description of the 
model is given ( for details see @ ). Section 3 deals with the implementation of 
the model. Section 4 describes the tuning setup and strategy. Results of tuning are 
summarized in Section 5. 
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2 (Modified) helix string model 



Assuming the string field can be modelled by a stream of soft gluons (Fig.Q]), direct 
hadrons acquire their transverse momenta by integration over soft gluon momenta 
along the field: 



A<J>, 

\p T \=2r\sin—\ (1) 

where pj is the transverse momentum of 
the hadron, r stands for the radius of the 
helix, and A<1> is the difference of the he- 
lix phase between the string break-ups 
which define the hadron. 
In the modified helix model, the helix 
winding is proportional to the energy 
density of the string: 




Figure 1 : Helix-like ordering of 
soft gluons at the end of partem 
cascade (left) and direct hadron 
transverse momentum (right). 



A® = SKAl = SE had (2) 

where K stands for the string tension [GeV/fm], A/ is the length of string piece, and 
S [rad/GeV] is a parameter. corresponds to the energy of the direct hadron in 
the string c.m.s. Parametrization (O implies strong correlation between energy and 
transverse momentum of the direct hadron (in the rest frame of the string). 

3 Implementation of the model 

The modified helix model is implemented in Pythia 6.421 |4] via private version 
of the fragmentation routine PYSTRF (6J. The following Pythia parameters and 
switches are used for steering: 

• MSTU(199)= 0/2 (standard fragmentation/modified helix fragmentation) 

• PARJ(102)= * ( helix radius r [GeV/c], replaces PARJ(21) ) 

• PARJ(103)= * ( helix radius variance [GeV/c] ) 

• PARJ(104)= * ( parameter S [rad/GeV] ) 

• PARJ(109)= 0.001 (azimuthal angle tolerance in the iterative search of the 
string break-up solution conform to Lund fragmentation rules and helix string 
parametrization) 
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4 Tuning setup and strategy 



The tuning is performed using the Rivet and Professor packages Q, using alter- 
natively Pythia pj ordered shower O, and Ariadne parton shower [8], on top of 
hadronic Z° decay generated by Pythia 6.421. 

The set of 6 simultaneously tuned parameters consists of 

• helix radius - PARJ(102) 

• helix winding - PARJ(104) 

• Lund parameter a - PARJ(41) 

• Lund parameter b - PARJ(42) 

• effective coupling constant Aqcd - PARJ(81) 

• parton shower cut-off - PARJ(82) 

(the later 2 parameters are replaced by PARA(l) and PARA(3) in the Ariadne tune). 

The setup of other Pythia parameters is based on the Professor tune of pj or- 
dered parton shower @, or on the DELPHI tuned setup 0, with small modifica- 
tions described in the Appendix. 

For simplicity, the helix radius variance is fixed (PARJ(103)=0.1 GeV/c). 

The tuning is performed on the set of inclusive charged particle spectra and 
event shape distributions measured by the DELPHI Collaboration Q. 

The input samples, generated with Rivet and private Pythia library, contain 
500k events each. The set of fitted optimal parameters is used to generate a vali- 
dation sample which is then compared to the reference data (i.e., the reported x 2 
difference refers to the validated result of the tune rather than interpolated esti- 
mate). 

5 Results and comparison with data 

Data distributions included in the fit are listed in Figf2] (for definitions see ll5l). 
The comparison of data with tuned model predictions is quantified by % 2 per bin 
measure, separately for each observable: 
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Figure 2: Comparison of tuned model predictions with DELPHI data, in % 2 per 
bin. Bar charts stand for modelling with Pythia pj ordered shower, points mark 
modelling with Ariadne parton shower. On average, helix string model (dashed 
chart and closed points) provides a better description of the DELPHI data. 
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With implementation of the helix string model, the string fragmentation effec- 
tively looses one degree of freedom, and the additional constraint translates into 
modified transverse momentum distribution. The change is visible in the low pj 
part of the transverse momentum distribution ( < 1 GeV/c ), and it is supported by 
the data (see plots on top of Figs. [3j Hit- 
Even more significant for the model validation is the improvement in the de- 
pendence of the average transverse momentum on the scaled momentum of the 
hadron (middle plots in Figs. OH]). This may well be the best evidence in favour of 
the existence of inherent correlations between transverse and longitudinal compo- 
nents of hadron momenta. 

The positive effect of the helix string modelling is seen also on a number of 
event shape variables (most notably on Major, minor; hemisphere broadening; dif- 
ferential jet rates). The overall % 2 for the whole set of distributions used in the 
tune is given in Table [T] and the numbers clearly indicate the preference of the 
data for the helix string model, when compared with the standard fragmentation of 
the Lund string. For completeness, the % 2 is quoted also for the set of identified 
particle rates, and the b-fragmentation function, even though these distributions 
are not included in the present tune. The tuned parameter values are quoted in the 
Appendix. 



Data set 


Pythia [ 10] 


helix + Pythia 


Ariadne 


helix + Ariadne 


inclusive spectra 
+ event shapes 
Nun = 619 


4075 


2453 


2453 


1489 


ident.part.rates 
+ b-fragmentation 

N bin = 47 


444 


669(*) 


614(*) 


586(*) 



Table 1: Sum (over all bins) of % 2 difference between data and models. 
The 'Pythia/Ariadne' labels distinguish between Pythia 6.421 pT-ordered parton 
shower, and Ariadne 4.12 parton shower. (*) distributions not included in the tune. 



Despite the overall improvement of the data description due to the helix string 
model, there are few distributions where a degradation is observed. There is some 
negative impact on the rapidity distributions (Figj5]>, and some degradation of the 
description of the scaled momentum is observed in the modelling based on Pythia 
pT-ordered parton shower (bottom plots in Fig©. The later is not confirmed in the 
study running Ariadne parton shower (FigJU), and there are some indications the 
effect may be related to the (untuned) charm quark fragmentation function. 
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Figure 3: Comparison of inclusive charged particle distribution measured by DEL- 
PHI and MC simulation based on Pythia pT-ordered shower, using standard frag- 
mentation ('Pythia 6.421') or helix string model ('helix+Pythia 6.421'). 
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Figure 4: Comparison of inclusive charged particle distribution measured by DEL- 
PHI and MC simulation based on Ariadne parton shower, using standard fragmen- 
tation ('Ariadne 4.12') or helix string model ('helix + Ariadne 4.12'). 
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Figure 5: Comparison of rapidity distribution ( w.r.t. sphericity axis ) as measured 
by DELPHI, and MC simulation based on Pythia pT-ordered shower (left) and 
Ariadne parton shower (right), for standard fragmentation (dotted line) and for 
helix string model (solid line). 



6 Conclusions 

A significant improvement in the description of hadronic Z° data is achieved with 
the help of modified string model which assumes the gluon field is ordered in a 
helix-like structure with pitch proportional to the energy density stored in the field. 

The improvement is achieved despite the minimal tuning setup which does not 
readjust the parameters of the heavy quark fragmentation. This means there is 
definitely room for further improvement of the tune. 

Tuned helix string parameters suggest a helix radius r ~ 0.4 GeV/c (rather well 
constrained) and a helix winding of S ~ 0.7 rad/GeV ( with large uncertainty ). 
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Appendix 

i/ Tuning of the model with Pythia p^-ordered parton shower 

The setup of the helix string model tune corresponds to the Professor flavour and 
p T tune with the exception of PARJ(47), fixed at 0.873 instead of 0.8. The 
change is done in order to facilitate comparison of the helix string tune with the 
updated Professor fragmentation tune ( iflOl ). 

Tuned parameter values and fit quality are given in Table [2] 



Parameter 


Tuned value 


PAPJ(41) 


0.084 ±0.12 


PAPJ(42) 


0.375 ±0.19 


PAPJ(81) 


0.237 ± 0.025 


PAPJ(82) 


0.65 ± 0.59 


PARJ(102) 


0.362 ± 0.032 


PARJ(104) 


0.510 ±0.23 


X 2 /( n bin-nd.o.f) 


4. 



Table 2: Results of the 6 parameter tune of the helix string model combined with 
Pythia 6.421 pT-ordered shower. 
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ii/ Tuning of the model with Ariadne parton shower 

The setup is similar to the one used in the DELPHI tune [5 ] with following modifi- 
cations: the detailed identified particle rate steering is not available, Bose-Einstein 
effect simulation is switched off ( MSTJ(51)=0 ), and radiative corrections are in- 
cluded ( MSTJ(107)=4 ). The tuning is performed twice, for the standard Pythia 
fragmentation, and for the helix string fragmentation. The tuned parameter values 
and fit quality are given in Tables |3j HI 



Parameter 


Tuned value 


PARA(l) 


0.240 ± 0.021 


PARA(3) 


0.537 ± 0.22 


PARJ(41) 


0.41 ±0.13 


PARJ(42) 


0.85 F 


PARJ(21) 


0.362 ± 0.032 


X 2 / ("bin - n d . .f) 


4. =2453/(619-5) 



Table 3: Results of the 5 parameter tune of Ariadne 4.12 with standard Pythia 
fragmentation. 



Parameter 


Tuned value 


PARA(l) 


0.264 


± 


0.024 


PARA(3) 


1.39 


± 


0.35 


PARJ(41) 


0.540 


± 


0.090 


PARJ(42) 


0.576 


± 


0.108 


PARJ(102) 


0.429 


± 


0.028 


PARJ(104) 


0.682 


± 


0.165 


X 2 / (n b in ~ n d .o.f) 


2.4 = 1489/(619-6) 



Table 4: Results of the 6 parameter tune of Ariadne 4.12 with helix string frag- 
mentation. 
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